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Abstract 

Modeling studies consistently demonstrate that the most violent winds in 
tornadic vortices occur in the lowest tens of meters above the surface. These 
velocities are unobservable by radar platforms due to line of sight consider- 
ations. In this work, a methodology is developed which utilizes parametric 
tangential velocity models derived from Doppler radar measurements, to- 
gether with a tangential momentum and mass continuity constraint, to esti- 
mate the radial and vertical velocities in a steady axisymmetric frame. The 
main result is that information from observations aloft can be extrapolated 
into the surface layer of the vortex. The impact of the amount of information 
available to the retrieval is demonstrated through some numerical tests with 
pseudo-data. 
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1. Introduction 

The strongest wind speeds in tornados are believed to occur a few tens of 
meters above the surface. Due to line of sight limitations, radar platforms are 
typically unable to measure this portion of the atmosphere. The relationship 
between the measurable flow aloft, and the unobservable (by radar) flow near 
the surface is complex (see for instance [2j, |7j, |9j and [TT] for different flow 
regimes). 

The reviews [10], [12] and [H] discuss the dynamics of different sections 
of a tornado. Snow ([E]) describes the change in magnitude of the different 
wind components both in the vertical and radial directions, which is based 
on simulations in fluid dynamics models and in the Tornado Vortex Chamber 
[1] at Purdue University. A tornado with a positive vertical velocity along 
the central axis is called a "single celled" vortex. The tangential velocity 
mean field increases as a function of height from ground level to a maximum, 
and then decreases again to the top of the vortex. Similarly, the tangential 
velocity increases as a function of the distance from the center of the vortex 
until it reaches a maximum, and then decreases to zero. This behavior can be 
captured with empirical parametric models, such as those discussed in [T5] . 
Models of this type have also been used in observational studies such as [16] 
to better understand measurements in the presence of noisy observations. 

In this paper, we estimate the three components of the wind velocity near 
ground level from observations aloft. The paper is divided into sections as 
follows. In Section [2| we review the basic considerations regarding observa- 
tions of atmospheric circulations by radar instruments and define the problem 
domain and relevant parameters of interest. Section [3] introduces a method 
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for estimating the vortex radial and vertical velocities, and Section ?? dis- 
cusses the mathematical issues related to this method. The mathematical 
issues include positive aspects of estimating flow fields with these dynamics, 
as well as situations in which the dynamics are insufficient to estimate the 
flow on the entire domain. Section |5] examines a few physical limitations of 
the approach. In Section |6| we perform an identical twin experimental test of 
the method for a tornado-like vortex. We generate pseudo-observations with 
an assumed tangential velocity model and random errors. Then we estimate 
the flow using the same tangential velocity model. This test is not meant to 
prove conclusively that the method will work with a real data set, but rather 
to show the theory in action. 

Remark 1.1. Many researchers in meteorology currently use variational 
techniques to estimate wind fields from radar velocity measurements. These 
techniques are powerful, and are especially useful for dealing with noisy mea- 
surements. They face the problems common to all optimal estimation tech- 
niques. Some of the difficulties are finding a unique global minimum and 
minimizer, and the tendency of least squares techniques to reduce the mag- 
nitude of smaller scale features. Further, a minimizer of a set of weakly 
enforced constraints may not satisfy any of the constraints particularly well. 
Boundary conditions for these types of methods are usually not chosen phys- 
ically, but rather are allowed to be retrieved with the rest of the variables. 
The authors are well acquainted with these techniques, and propose the tech- 
niques in this paper as a first step toward remedying some of these difficulties. 
Most variational techniques utilize some sort of descent based minimization 
procedure, and the solutions provided by the method in this paper could be 
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used as the "first guess" which is required of all iterative schemes. 
2. Background 

Assume that two radar instruments measure a given volume of air simul- 
taneously. The two horizontal components of the velocity can be recovered 
if the radar beams are approximately horizontal. In this case, the mea- 
surements contain very little information about the vertical component of 
velocity, i.e. are orthogonal to the vertically pointing basis vector. Take the 
flow to be in cylindrical coordinates, with the axis of the coordinate system 
aligned along the vertical axis of the vortex. Thus the recovered components 
are the tangential and radial components of the swirling flow. 

For the remainder of the this work, assume two sets of wind measure- 
ments, which have been converted to radial and tangential velocities for the 
vortex of interest, and averaged azimuthally to create an axisymmetric mean 
pair of velocities. The spatial domain includes the vertical axis and the sur- 
face and measurements which are representable by a parametric model. A 
family of parametric models for the tangential velocity is chosen which best 
approximate the qualitative features of the given data, then a particular pa- 
rameter set is selected so that the tangential velocity model is optimal (in 
some sense). This is done in advance of seeking to estimate u and w. 

In the next section, the estimation of radial and vertical velocities in a 
layer near the surface, where the velocities are not observable, is considered. 
The problem is posed on the domain Q, which is illustrated in Figure [l| The 
domain is decomposed into an observable region Qo and an unobservable 
region Q^, separated by a horizontal line z = h. This line is referred to 
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Figure 1: Schematic of Problem Domain. The outer radial boundary is dashed to represent 
the unknown boundary condition. 

as the minimum observable height (MOH) line. The domain on which we 
interested in retrieving the flow is referred to as the surface layer, which is 
the portion of the domain between the height z = and z = kg, where we 
will refer to hs as the surface layer height. The parameter hg is chosen for the 
application of interest. For example, if we are interested in surface damage, it 
might suffice to only examine the flow in the layer with hg = 1 meter, whereas 
structural engineers might be interested in multistory buildings, and would 
necessarily use a larger value for this parameter. 

3. Estimating u and w 

Assume that the vortex is approximately steady and axisymmetric, and 
that f (r, z) captures the essence of the tangential velocity present in the 
observations. Consider the steady, axisymmetric Navier-Stokes equations of 

5 



motion, given by 
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where u, v and w represent the three components of the velocity vector in 
cyhndrical coordinates, p the density, p the pressure, and u the fluid viscosity. 
Further, assume that the fluid is incompressible, and so mass conservation 
takes the form 
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If v{r,z) were a component of a true solution to ([T])-([3]), then this system 
would still have a unique solution if v exactly modeled error free data, and 
if these equations exactly hold for real atmospheric vortices. Realistically, 
observational and model errors lead us to conclude that enforcing only a 
subset of these dynamics may help to avoid an overdetermined problem. 

Introduce the vertical vorticity Cl'"?^) = r^^ar^ radial vorticity 

ri{r,z) = — If. With these substitutions, ^ can be rewritten as 



z)u{r, z) — rj{r, z)w{r, z) = u 
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which is an algebraic relation between u and w, once v has been selected. 
Next, introduce a streamfunction \E', defined by 
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in cylindrical coordinates, so that \1/ satisfies ^ automatically. The tangen- 
tial momentum equation ^ becomes 



C{r, z) — {r,z) + r]{r,z) — {r, z) = z/r 



dr 



oz 



r, z) — —[r, z) 



(6) 



This is a hyperbolic boundary value problem on Vlh- The boundary condi- 
tions at the surface and vertical axis should yield vortical flows similar to 
actual atmospheric vortices. By choosing zero Dirichlet boundary conditions 
for \1/ on the lower and axial boundaries, mass is conserved. The bound- 
ary condition for u along z = h is provided by the measurements, and the 
boundary condition for w is taken to be the result of solving ([s]) for w and 
substituting in the condition for u. Once w(r, z) is known along the MOH 
line, \l/ is recovered using 



^(r, h) 



sw{s, h)ds. 
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The outer radial boundary is left unconstrained for the moment. 

Equation ^ is quasihnear with associated characteristic equations [6j: 
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where t is the characteristic variable for the position along the characteris- 
tic curve given by {r{t), z{t)). To seek solutions these ordinary differential 
equations must be supplemented with initial conditions. Let s denote the 
characteristic variable which distinguishes between different characteristic 
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curves, by parameterizing the initial values for r, z, and and define 



r(0,s) = s (11) 
z{<d,s) = h (12) 

^(0,s) = - f sw{s,h)ds (13) 
Jo 

This choice of initial conditions means that the equations are initialized with 
values on the upper boundary of Qh and allow the dynamics to propagate 
the information contained on them down into the domain. 

Remark 3.1. Assuming that v G C''{Qh) for some k > 2, classical results 
from the theory of ordinary differential equations (for example, those in [5]) 
provide existence and uniqueness of solutions to these initial value problems, 
and smoothness with respect to the initial conditions. This implies that if a 
point (r, z) lies on a characteristic curve that intersects the upper boundary 



of ilh, there is a classical solution \1/ defined at {r,z) that satisfies (10) and 



(13). 



Remark 3.2. The fluid viscosity u is an important physical constant for the 
purposes of time dependent model simulation. Since the flow is stationary, 
u has a small impact on the results with this method. Where it makes 
calculations simpler, u will be set to 0. 

In order to simplify the discussion, we introduce the following notation. 

Definition 1. For a point {ro,Zo) € Qh, define 

(1) c(-,ro,Zo) : M — 7- the solution mapping of the dynamical system 
with initial condition (ro,Zo). 



(2) C{ro,Zo) = c{M.,ro, Zo) : the set of all points {r,z) which can be at- 
tained by integrating ([8])-(|9| (either forwards or backwards) starting 
from {ro,Zo). 

(3) Kh = {{ro.Zo) G nu\{r,h) ^ Cijo.Zo) Vr G [0,i?]}. 

C(r, z) is referred to as the characteristic curve containing {r, z). The set 
is referred to as the information void for the problem, because the dynamics 
do not carry information from aloft to these points. 

4. Surface Layer Wind Velocities 

Assume that v{r,z) = (f){r)tp{z), where 

Assumption 4.1. (i) 4> andip both are k times continuously differentiate 
{k > 2). 

(2) (no-slip condition) 0(0) = ^'(0) = 0. 

(3) (j)>0 on {0,R) anc?^ > on {0,H). 

U) fri^o) + i<t>{ro) =0 and'-'-^j^Oforrj^ r,. 

(5) ^^{z,) = and'^y^Oforz^z,. 

This assumption allows a more thorough analysis, and [15] has demon- 
strated the utility of such models for data analysis. A schematic streamfunc- 
tion of a vortex embodied in these assumptions is shown in Figure [2] 



The following result says that Assumption 4.1 always yields a nontrivial 
surface layer in which we can retrieve the flow: 
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Figure 2: Schematic of Characteristic Curves under Assumption |4.1[ where the relative 
maximum of F occurs at (ro,Zo). 



Theorem 4.2. // Assumption 4-1 holds, then there is an ho such that if 
z<ho, C(r,z)n(fi\fi/,) ^0. 



The next four resuhs form the basis of the proof of Theorem 4.2 



Lemma 4.3. The characteristic curves C{r,z) are the level curves of T, 
where T = rv is the circulation on circles parallel to the horizontal plane, 
centered on the vertical axis. 

Proof. When 77 7^ 0, we can write the solution curves as {r,z{r)) by consid- 
ering 

ar rj 

and when ( ^ a.s {r{z), z) from 

dr rj 
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Note that 

C _ {rv)r _ 

which imphes that the characteristic curves are everywhere tangent to the 
level curves of F. Hence, viewed in the plane, these two collections of curves 
are the same. □ 

Remark 4.4. This result specifies the characteristic curves in terms of our 
tangential velocity model, which is estimated a priori utilizing a least squares 
(or some other) data mismatch criterion. It also gives a criteria by which 
to avoid the solution ^ being multiply defined, which can occur when using 
method of characteristics. To avoid this behavior, choose v to be appropri- 
ately smooth. 

The next result states that when the maximum tangential velocity is in the 
observable region, then the flow is retrievable over all of Q using character- 
istics. 

Lemma 4.5. If h < Zo, then Kh — 0. 

Proof. There are two cases. For r < Vg, r] < and C > 0. Hence if C(r, z) is 
traversed in the positive t direction, the curve must eventually cross z — h, 
since C(r, z) cannot intersect the vertical axis. For r > ro, rj < and C < 0. 
Since C (r, z) cannot intersect the horizontal axis, there must be a Zi such 

that {to, Zi) G C(r, z), and now apply the argument from the first case, using 
(ro, zi) as our initial point. Hence, for any (r, z) e flh, C{r, z)n{fl\ flh) 7^ 0, 
which implies Kfi — $. □ 
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Lemma 4.6. Suppose that < zi < H satisfies C{ro,zi) C Q°. Then 
C{ro,zi) is a closed curve. 

Proof. Let zi > Zo- Since F has only a single relative maximum, |^ > 
for r < To and < for r > Tq, and similarly for the vertical gradient 
of F. Consider the characteristic curve which passes through (ro,zi), and 
first traverse in the positive t direction. Since Vz{ro,zi) < 0, > 0, and so 
the characteristic curve moves to the right. For r > Tq and z > Zo, rj > Q 
and C < 0, and so the characteristic curve moves to the right and down. 
Since C(ro,-2i) C there must be an ri with < ri < R such that 
(ri, Zo) G C(ro, zi)^ else C(ro, zi) would cross the line r = R. Similarly, since 
for r > To and z < Zo, rj < and C < 0, there must he a < Z2 < zq such 
that (ro, ^2) G C{ro,zi). Otherwise C{ro,zi) would intersect the lower axis 



z = Q, which would contradict Corollary 5.1 Thus C{ro,Zi) intersects the 
line r = To at {ro,Z2). By traversing C{ro,Zi) in the negative t direction 
starting from (ro,-2i), and using similar arguments, there is a < z^ < Zo 
such that (ro, ,23) G C{ro,Zi). 

Suppose Z2 < Z3. Then there is a z* with Z2 < z* < Z3, and since F^ > 0, 
we must have that 

T{ro,Z2) <T{ro,z*) <T{ro,Zs). (16) 

But this is a contradiction, since F(ro, Z2) = F(ro, ^3). 

Let t2 such that c(t2,ro,Zi) = {ro, Z2) and ^3 such that c(— 13, To, -Zi) = 
(ro, Z3). Then c(t2 + ^3, Zi) = {ro, zi) and C{ro, Zi) is closed. 

If zi < Zo, a. symmetric argument shows that C{ro, zi) is closed. □ 
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Lemma 4.7. Suppose that Zo < h. Then one and only one of the following 
statements holds: 



(1) C{ro, h) is a closed curve, and Kh is the interior of the region enclosed 
by C{ro, h). 

(2) C{ro, h) intersects the outer radial boundary at [R, zi) and [R, Z2), and 
Kh is the interior of the region enclosed by C{ro,h) and the segment 
{{R, z) : zi < z < Z2} 

Proof. First, if C{ro,h) is not a closed curve, then if C{ro,h) is traversed 
in the negative t direction, it must cross the hne z = Zo, and then the hne 
r = To, because C{ro,h) cannot intersect the axes. This imphes that there 
is a t such that c{—t,ro,h) = {r*,z*) with r* > To and z* < Zo- If C{ro,h) 
were to cross the hne z = Zq again, then the signs of the vorticities would 
force C(ro, h) to intersect r = Tq, and at the point (ro, h) by the argument in 



Proposition 4.6 Similarly, if C(ro, h) is traversed in the positive t direction, 
C(ro, h) cannot cross the line z = Zo, or else C(ro, h) would be a closed curve. 
Thus, either C(ro, h) is a closed curve, or C(ro, h) intersects the outer radial 
boundary at two distinct points {R,Zi) and {R,Z2), where Zi < Zq < Z2. In 
either case, denote the set enclosed by C(ro, K) (and possibly {i?} x [^1,^2]) 
byifo. 

If (r, G \ -K^o) proceed as before by traversing C(r, either in the 
positive (r < To or z < Zq) or negative (r > Tq and z > z,^ t direction. Note 
that c(t,r, z) ^ Ko for all t G M because dKo = C{ro,h) (possibly plus the 
outer boundary), and characteristic curves may not intersect. Since c{t,r,z) 
also cannot intersect the axes, there must be a t such that c(t, r, z) E Qh)- 
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Thus K C Ko. 

If (ri,zi) G Ko, then r(ri,zi) > T{r,z) for all (r, z) E Q \ Qh- Thus 

C(ri,2:i) n = and so (ri,2;i) G -ft'. Hence Ko C -ft', and K = 

Ko. □ 



Proof of Theorem 4-2 



Proof. If Zo > h, simply take ho = hhj Lemma 4.5 Suppose Zo < h.d Since 
C{ro, h) C Q is either closed or intersects the outer radial boundary, it is also 
compact. Hence the map (r, z) ^ z has a minimizer at some point ho. Thus, 



if z < ho, (r, z) ^ Kh, and so C{r, z)n{Q\ Qh) 7^ by Proposition 4.7 □ 



Remark 4.8. The height z = ho can be referred to as the minimum un- 
reachable height, since for values of z < ho, the solution is reachable via 
characteristic curves. Under Assumption 4.1[ the proof of Theorem 4.2 



im- 



plies that ho is the smallest solution of ipi^z) = ipih). 

Remark 4.9. A similar result holds for the map (r, z) t— )■ r, implying the 
existence of a "minimum unreachable radius", though this is not directly 
relevant to the problem initially posed. 

5. Model Limitations 

5.1. Boundary Conditions 



The following corollary follows immediately from Assumption 4.1 and Lemma 4.3 



Proposition 5.1. If Assumption 4.I holds, then no characteristic curve may 



intersect the lines r = and z = 0. 
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Proof. Note that r(0, z) = r(r, 0) = 0. Since F > on the interior of Qh, 
no level curve of F intersections the boundaries, and hence no characteristic 
curves intersect the axes. □ 

Remark 5.2. Proposition |5. 1 1 implies that the choice of boundary conditions 
along the surface and the vertical axis do not affect the flow on the interior 
of the domain, so long as v vanishes on these axes. This is a consequence 
of the choice of dynamical constraints, and removes a physical degree of 
freedom from the problem, since in real vortices, surface roughness effects 
can propagate into the domain. The literature contains multiple discussions 
(e.g. [Z], [in]) of what boundary conditions are most realistic, and generate 
physically realistic vortices. It is intuitively clear that the radial and vertical 
velocities will depend on the their behavior at the surface and along the 
central vertical axis, but this is not captured by the dynamics we are choosing 
to constrain the solution. 

5.2. Multiple MOH Intersections 

Another difficulty is the possibility of characteristic curves intersecting 
the MOH line multiple times. In this case, the boundary data on the MOH 
line may not be compatible with the dynamics. For real data, this will al- 
most certainly not be the case due to noise and the error introduced by the 
tangential model v. This situation is reminiscent of the data assimilation 
problem that is usually tackled using least squares minimization of an ob- 
jective functional that penalizes disagreement between model prediction and 
observation relative to the uncertainty present in each. More information 
about this topic is found in [Sj. This problem will be addressed in a future 
work. 
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5.3. Velocities above the Minimum Unreachable Height 

The results in Section |4] point to potential difficulties when kg > ho, 



the minimum unreachable height guaranteed by Theorem 4^ Clearly, there 
are portions of this set that are reachable by characteristics, namely those 
characteristic curves that pass through to the surface layer below ho- The 
rest of Qh is precisely K^, which we have called the information void. 

6. Numerical Experiments 

In this section, a simple test of the theory developed in Section |4] is 
demonstrated. This experiment is an identical twin, since the same functional 
form is used to generate the observations as the one used to select v and 
estimate u and w. Initial tests showed that dependence on the viscosity u 
was small. With this in mind, assume z/ = 0, which simplifies the numerics 
from solution of a linear ordinary differential equation for $ to solving the 
equation $(r, z) = $(ro, This equation can be approximately solved to 
any specified degree of accuracy using a simple bisection method. 

6.1. Generation of Pseudo-observations 

As a first experiment, a collection of pseudo-observations is generated 
that emulates a single time of model output from Davies- Jones' axisymmet- 
ric model, described in [1]. At the time of interest, the tangential velocity 
near the surface exhibits a single maximum. The radial velocity is negative 
beneath this maximum, which is typical of a swirling flow (|13]) with a no- 
slip condition on the tangential velocity, and represents air being drawn into 
the vortex. Finally, the vertical velocity is relatively large and positive along 
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the axis adjacent to the tangential maximum, which is also typical of these 
types of flows. The tangential velocity v{r, z) is modeled using a product of 
functions of the form 



ri — 1 



(j)^w{x,n,Xc) = . — -. (17) 

[n — IjXc + 

The function cj) has a smooth maximum at (xc, 1), and increases approxi- 
mately linearly on (0, Xc), and decays like x^~^ as x co. Assume 

v{r,z) = Vc4>ww{r,nr,rc)(l)ww{z,nz,Zc). 



This function satisfies Assumption 4.1, and so all of the theory in Section ?? 



is valid for this choice of model. The velocity pseudo-obs used are depicted in 
Figure [3| with the tangential velocity depicted as contours, and the radial and 
vertical velocities depicted as a single vector in the r — z plane. To simulate 
the effects of measurement error, a set S of independent realizations of a 
normal distribution with standard deviation equivalent to 1 ms~^ was added 
to the tangential velocity values on a discrete spatial grid. Three separate 
experiments were run to simulate random errors with deviation as large as 3 



ms ^. 



6.2. Impact of MOH on Surface Layer Thickness 



As a demonstration of Theorem 4.2, the streamfunction for a fixed initial 



condition was computed for MOH values of 1.5, 2.5 and 3.5. The resulting 
surface layer streamfunction is plotted in Figure |4} Note that as more of the 
vortex is observable, more is retrieved below the MOH line. Also note that 
even in the case with the least information (h = 3.5), there is a retrieved 
surface layer of nontrivial thickness. 
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Velocity Pseudo-observations 




0.5 1 1.5 2 2.5 3 3.5 4 

r 



Figure 3: The "truth" velocity fields, where the contours are of tangential velocity with 
increment 0.1 from to 1, and the vectors indicate the radial and vertical velocities, with 
a reference vector of length 1 shown in the bottom right corner. Note the strength of the 
infiow near the surface, and of the updraft near the vertical axis. 

7. Discussion 

A methodology was introduced for extrapolating observations of wind ve- 
locities downward toward the surface. For the dynamics chosen to constrain 
the flow, the information contained in observations aloft propagates along 
curves that coincide with the level curves of F = r^;, which is estimated from 
observations in advance of the problem discussed here, and hence is known 
a priori. With more assumptions about our tangential model, the location 
and size of the information void Kh are exactly known for a specific value of 
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h = 3.5 h = 2.5 h=1.5 




Figure 4: Comparison of Retrieved Surface Layers for Three Different MOH Values 

h. An important result is that there is always a nontrivial height /lo, below 
which, everything can be retrieved using the characteristic framework. 

As a first test of the method, a set of pseudo-observations was created 
that approximate the behavior of simulated tornadoes. The fiow was then 
estimated assuming different amounts of knowledge, embodied in the param- 
eter h, and different levels of random error. These experiments yield mixed 
results, because it always occurs that the true value of , , or ItTlmax falls 
inside the spread of results, but not always with the correct frequency (rela- 
tive to agreement with the tangential velocity observations). It is clear from 
the theory that the functional form chosen for f (r, z) has strong impacts on 
what can be retrieved, as well as the quality of what is retrieved. 

The authors assert that knowledge of the shortcomings of a particular 
method is valuable information; hence, the analysis of and focus on informa- 
tion voids. In a more standard variational technique, these voids would not 
appear, because the various smoothness terms would ensure that a smooth 
solution is defined everywhere. However, the solution in those regions is no 
more physically relevant than any other solution, since it is completely de- 
termined by the terms that are introduced for numerical stability. A natural 
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next step is the inclusion of another dynamical constraint, such as the bal- 
ance equation for azimuthal vorticity. This extra constraint could yield a 
method to estimate a meaningful solution in the information voids, without 
relying on unphysical smoothing. 
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